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We have investigated temperature trends of the microscopic structure of the SPC/E and TIP4P-Ew water models 
in terms of the pair distribution functions, coordination numbers, the average number of hydrogen bonds, the 
distribution of bonding states of a single molecule as well as the angular distribution of molecules by using 
the constant pressure molecular dynamics simulations. The evolution of the structure is put in correspondence 
with the dependence of water density on high temperatures down to the region of temperatures where the 
system becomes supercooled. It is shown that the fraction of molecules with three and four bonds determine 
the maximum density for both models. Moreover, the temperature dependence of the dielectric constant is 
obtained and analyzed. 
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1. Introduction 

It is our honor and pleasure to contribute to this special issue and to dedicate this paper to Prof. 
Dr. Douglas Henderson on the occasion of his 80* birthday. Almost four decades ago, Barker and Hen¬ 
derson published their seminal paper “What is “hquid”? Understanding the states of matter”, that com¬ 
prehensively describes basic theoretical methods to explain the equilibrium properties of simple liquids 
Iin . These methods have made substantial progress since that time. In particular, computer simulations 
techniques have become essentially powerful and common tools in the theory of liquids, mixtures and so¬ 
lutions. Computer simulations permitted to reach enormous progress in understanding the liquids with 
hydrogen bonds, more specifically, water and aqueous solutions, for needs of physical chemistry and for 
several inter-disciplinary areas involving complex biological molecules and membranes. In particular, 
Henderson with co-workers was actively interested in the properties and behavior of inhomogeneous 
aqueous solution, see e.g. Sfll. 

The studies of water have long history that has been documented in an enormous number of publica¬ 
tions. Of particular interest are the peculiar thermodynamic and dynamic properties of water manifested 
in a set of anomalies. In spite of formal equivalence of the relation between the microscopic structure and 
thermodynamic properties discussed in great detail by Barker and Henderson for simple homogeneous 
liquids |l|], in the case of water it is inevitable to involve the notion of hydrogen bonding. Many works 
have dealt with the exploration of hydrogen bonds network structure, its dynamics and life-time of bonds 
at an early stage of computer simulations of aqueous systems, see e.g. a few original papers @,01 as well 

* E-mail: pizio@unam.mx 


© E. Galicia-Andres, H. Dominguez, 0. PizioE-mail: pizio@unam.mx, 2015 


13603-1 





E. Galicia-Andres, H. Dominguez, 0. Pizio 


as later reviews , that contain an extensive list of references on the subject. Moreover, it is known 
that unique properties of water specifically manifest themselves in the region of states below the freezing 
point corresponding to the supercooled liquid, see e.g. a recent review lldl and references therein. 

It is seducing to relate at least some of water anomalies with the transformations occuring in the 
hydrogen bonds cooperative structure. Therefore, in this short report our focus is in the description of 
the evolution of the average number of H-bonds per water molecule and in the distribution of the num¬ 
bers of H-bonds per molecule upon temperature changes. Moreover, we would like to relate the trends of 
the behavior of these properties with the density anomaly. While the properties in question describe the 
formation of “agglomerates” of bonded molecules at relatively short-range distances between them, we 
would like to get additional insights into the simultaneously occuring changes of the long-range correla¬ 
tions between water molecules in terms of the dielectric constant. 

One of the most studied anomalies of water is its nonmonotonous temperature dependence of density, 
p{T), see e.g. a review providing a nice historical insight j^. Several non-polarizable water models with 
different number of interaction sites are capable of reproducing the presence of maximum density of 
water with different precision. 

Namely, the temperature of maximum density (TMD) at ambient pressure for most frequently used 
models is observed at 2S3 K (TIP4P), at 278 K (TIP4P/200S), at 274.1S K (TIP4P-EW), at 241 K (SPC/E) and at 
182 K (TIP3P). These values for temperature are taken from the table 2 of the review by Vega and Abascal 
illll . The experimental results for these and other properties used by the authors are taken from jlill . 
The experiment yields 277 K for the TMD of real water. On the other hand, the water density at TMD is 
0.988 (TIP4P), 0.993 (TIP4P/2005), 0.994 (SPC/E), 0.98 (TIP3P), while the experimental value is 0.997. Some 
deviation from these values reported in different works result from technical details of each simulation. 

The origin of density anomaly is commonly attributed to the evolution of hydrogen bond network 
with the formation of an open structure as the temperature is lowered or, on the other hand, manifests 
itself as the reflection of a possible liquid-liquid phase transition observed by simulations in some water 
models Ql. The latter and other scenarios for peculiar properties of supercooled water were documented 
in{]3]. 

To summarize this brief introduction, our intention in this work is to to present our own data of 
molecular dynamics simulations of the temperature dependence of water density, p{T), in an ample in¬ 
terval of temperature (from 370 K down to 170 K) at ambient pressure by using the SPC/E model. Next, we 
would like to compare the obtained data with the results previously reported by other authors. However, 
the principal focus is to follow changes of water structure in terms of different descriptors involving both 
the hydrogen bonds and the p{T) dependence in order to establish a clear relationship between them. 
Nevertheless, for the sake of comparison of our observations for the SPC/E, and to make the conclusions 
more sound, similar results for the TIP4P-Ew water model lllill are presented as well. Our interest in 
the SPC/E model is motivated by its wide use for different purposes, in particular for the description of 
thermod)mamic properties of mixtures. Extending the previous research in mixtures, having in mind a 
more ample project presently being developed in this laboratory, our intention is to use the SPC/E model 
as a starting point in exploring water-alcohol and water-DMSO solutions by using molecular dynamics 
simulations. 


2. Simulation details 

All our simulations of water models were performed in the isothermal-isobaric ensemble at ambient 
pressure. The initial configuration was constructed with 2000 water molecules placed in a cubic array in 
a simulation box. Then, the system was simulated in an ample range of temperatures, from a high tem¬ 
perature 370 K down to 180 K. The DL-POLY Classic package was employed (l^. We used the Berendsen 
thermostat and barostat, the running timestep was set to 0.002 ps. Commonly, the periodic boundary con¬ 
ditions were used. The short range interactions were cut-off at 11 A, whereas the long-range interactions 
were handled by the Ewald method with the precision 10“^. After equilibration, several sets of simulation 
runs, each for 6-8 ns, with restart from the previous configuration, were performed to obtain averages 
for data analysis. The overall length of simulation was dictated by the stability of internal energy, density 
and the dielectric constant, all being the functions of time. 
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3. Results and discussion 

The temperature dependence of the density of water and the stability of several ice phases in the 
framework of the SPC/E model were studied in several works, see e.g. IlsU^ . In figure[T](a) we show the 
results obtained by different authors for the SPC/E water and the results of the present work. The liquid 
and glass branches, denominated in figurejl](a) as Baez et al., were constructed using the points reported 
in table II and in table III of reference (l8|]. The data follow from the simulations of 360 molecules. It 
was assumed that the cut-off distance is at 9 A and the effect of interactions at larger distances was 
neglected. The density maximum (p = 1.0265 g/cm^) was obtained at 235 K. In the later work from that 
laboratory, the density values for disordered hexagonal ice Uh) in the interval between 170 K and 220 K, 
were reported, see reference These data are also shown in figure[T](a). However, it is important to 
mention that the results refer to the modified SPC/E model, i.e., by multiplying the interaction potential 
by a switching function between 6 A and 9 A to make the forces and the potential continuous at a cut-off 
distance. The modification of the potential and the application of Ewald summation technique to include 
the long-range interactions also resulted in the changes of the estimates for TMD, in the improved 

values are p = 1.003 g/cm^ and T =; 260 K. The algorithmic peculiarities of these series of works include a 
profound exploration of the effects of Ewald summation on the temperature dependence of density and 
other properties of the system as well as calculations of the Gibbs free energy of distinct phases which 
requires thermodynamic integration. The line describing the density of stable /h ice obtained by Gay et 
al. is reproduced from fisll . It was mentioned therein that the stabihty of this specific crystalline phase 
crucially depends on the number of molecules in the simulation cell. It is worth mentioning that the 
authors performed simulations in isostress, constant temperature ensemble and estimated the melting 
point by analyzing the time evolution of the properties of the system upon heating of the sohd phase. 

Quite recently, Vega et al. (i^l evaluated the melting temperature of the most popular non-polarizable 
water models. Eor the SPC/E model, the melting temperature value of ice 1^ and the densities of liquid wa¬ 
ter and ice are shown in figure[T](a). Our data for the liquid density of the model starting from a rather 
high temperature down to the region of supercooled liquid are compared with the previously published 
results by Bryk and Haymet 11^ a nd with very recent results by Fennell et al. (2l|]. A small difference 
between the curve obtained in 11711 and the present results in the entire interval of temperatures is pre¬ 
sumably caused by technical details, namely by the distinct cut-off value. On the other hand, it is worth 




Figure 1. (Color online) Panel a: Temperature dependence of the SPC/E model water density in the liquid 
and glassy state as weU as of /j, ice at a constant ambient pressure. The results of different authors are 
marked in the figure. Panel b: A comparison of the density dependence of the SPC/E and TIP4P t5rpe 
models. Circles show our results for the TIP4P-Ew model. 
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mentioning that Fennell et al. jill] slightly modified the original SPC/E model by fitting the dielectric con¬ 
stant value at room temperature to the experimental result. Our data locate the TMD at 245 K, i.e., in 
between the commonly used value 241 K, see e.g. (nl[^, and the result for the optimized model by Fen- 
nell et al. 247 K. As for the value of density at TMD, our result is quite close to the data reported in flTII . 

In general terms, the water density values in relation to temperature grow in the interval that starts 
at high temperatures and is delimited below by the TMD. Next, if the temperature decreases further, the 
density decreases until it reaches a minimum value, Tnun. our data locate the minimum at T = 190 K. The 
minimum value of density was observed for a set of water models, e.g., for ST2 @] and TIP4P-200S (H). In 
the case of the SPC/E model, the melting temperature Tm is located in the interval between the TMD and 
Tmin. and the region between the melting point and Tnun actually corresponds to the supercooled states. 
Einally, if the temperature decreases below the Tnun. the density grows again as in a normal fluid, but this 
growth is not as rapid as in the interval of above the TMD. 

In our opinion, a plausible explanation of the presence of the density minimum that involves struc¬ 
tural changes in water was given in @|. It refers to the behavior of a simple liquid that shrinks upon 
cooling due to the increasing effects of attraction between particles with a decreasing temperature. How¬ 
ever, in a certain temperature interval, this normal behavior alters due to augmenting effects of hydro¬ 
gen bonding between molecules and the resulting formation of “complexes” or “agglomerates”. Coop¬ 
erative behavior of these entities result in the growth of an expanded, open structure. The effects of 
directional hydrogen bonding overcome the effects of an attractive interaction. Consequently, the fluid 
density decreases in a certain interval of temperatures, evidently specific for each model. However, the 
directional interactions are characterized by saturation. Therefore, upon further cooling, when changes 
of the structural units become weaker, compared with the previous regime (as documented for example 
by the changes of the average number of hydrogen bonds per molecule), then the fluid density grows 
again similar to the temperature interval above the TMD. 

However, the above discussion is incomplete and some thermodynamic arguments must be taken 
into account. The liquid is the stable equilibrium phase for temperatures above the melting temperature 
and it is unstable below a stability temperature, whereas in between, the liquid is metastable with respect 
to crystal ice (^. In general, metastable liquids and specifically water, can survive even at temperatures 
well below the melting point until homogeneous nucleation takes place and water converts into ice jiill . 
Therefore, a certain part of the curve p(T) presented in figure[T](a) for the SPC/E model (below the melt¬ 
ing point) characterized by a descending density can correspond to metastable, “long-fife” states (on tfie 
simulation time scale). Similar comments concern the behavior of the TIP4P-200S model, see figure 4 of 
reference and the ST2 model 

Due to finite-size effects and consequent details in implementation of Ewald sums as well as due to 
finite sampling time, a possible spontaneous crystallization, however, is not observed in many works 
for a set of water models, as well as in our present simulations. Debates concerning this problem just 
started very recently, see e.g. j^lUll • These authors provide a comprehensive discussion of various issues 
relevant to the problem using their own data coming from the Monte Carlo and molecular dynamics 
simulations and analyzing the previously published results for different water models from the literature. 
Some other papers from the special issue of the Journal of Chemical Physics devoted to confined and 
interfacial water published in November 2014 can be useful for the reader as well. 

Our final comments concerning a set of curves shown in figure [T] (a) refer to the glassy water in 
the framework of the SPC/E model. Actually, water, if considered below the homogeneous nucleation 
temperature, should be in a glassy form, see e.g., a rather comprehensive discussion of the issue in j^. 
In the case of the SPC/E model, the glass transition temperature was located around 177 K (l3l. At this 
temperature, the dependence of internal energy on temperature shows a kink (we reproduced the results 
of that work by performing long simulations with 2000 molecules and located the kink around T =; 185 K). 
Thus, the augmenting density region at low temperatures refers in part to the amorphous states that 
cannot be distinguished in terms of the structural characteristics such as the pair distribution functions 
coming from the present simulations. 

In figure [T](b), we present the temperature dependence of density of the SPC/E model and compare 
it with our results obtained for the TIP4P-Ew model performed in the framework of the above described 
methodology. In addition, the results by H. Pi et al. jH] for the liquid and ice 1^ phases are given as 
well. The TIP4P-Ew model reproduces the experimental result for the TMD and an overall dependence of 
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density at ambient pressure much better compared to the SPC/E. All the models involved in this figure 
(SPC/E, TIP4P-EW, T1P4P-2005) predict the existence of the density minimum and augmenting density 
at lower temperatures. However, it is worth mentioning that the density of ice 7h substantially differs 
from the density of supercooled water (and of its glassy form) for the SPC/E model [see figure [T] (a)], 
whereas in the case of T1P4P-2005 model, the licfuid phase density and that of ice are close to each other 
as it follows from the simulations by H. Pi et al. [2^, see figure [l](b). Similar behavior concerning the 
density dependences has been documented recently by J. Alejandre et al. (3 for fho model from the 
T1P4P family developed by fitting the dielectric constant, the maximum density and the equation of state 
at low pressure. Unfortunately, we are unaware of the simulation results for the 1^ ice phase of TlP4P-Ew 
model at room pressure to make our discussion more general at this point. However, it seems that the 
differences between liquid and ice phase densities in the region of low temperatures are very sensitive 
to the minor changes of parameters describing different models. 

The microscopic structure of a set of water models coming out from MD simulations with respect to 
the experimental diffraction data at room temperature and ambient pressure was discussed in great de¬ 
tail by Pusztai et al. by using the protocol combining the experimental total scattering structure 
factor and partial radial distribution functions from simulations. It was shown that the TtP4P-200S model 
exhibits the best consistency with the experimental structure compared to other commonly used mod¬ 
els. Nevertheless, the SPC/E, ST2 and TtP4P results are not seriously inconsistent at ambient conditions. 
Unfortunately, such an analysis of the structure of water has not been performed for other temperatures 
and pressures due to experimental difficulties. 

A few examples of the radial distribution functions, gtj, where i,j stand for 0 and H, obtained in 
our NPT simulations of the SPC/E model are presented in two panels of figure |2] It can be seen that 
the height of the first and the subsequent maxima of both functions, gooir) and gou. increases with a 
decreasing temperature in the entire interval of temperatures studied. These trends are well pronounced 
upon cooling in the range starting from a high temperature down to around 190 K, further. At even lower 
temperatures, the overall shape of gijir] is similar to that observed at 190 K. However, changes of the 
structure are much less pronounced if one moves to 180 K and 170 K. We omit these results for the sake of 
brevity. The functions gooir) and gonC^) are given mostly as an illustration, because their characteristics 
are necessary to introduce the concept of hydrogen bonds between molecules. 

In terms of the coordination number, noo. it is appreciable that the cooling causes more pronounced 
shells around the chosen molecule, cf. figure(a). While af high femperatures even fhe first coordina- 



Figure 2. (Color online) Oxygen-oxygen (panel a) and oxygen-hydrogen (panel b) pair distribution func¬ 
tions of the SPC/E model at different temperatures. 
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Figure 3. (Color online) Oxygen-oxygen coordination number of the SPC/E model at different temper¬ 
atures (panel a). The derivative of the oxygen-oxygen coordination number at different temperatures 
(panel b). 


tion shell is not very well pronounced, it is seen perfectly well at low temperatures, and the coordination 
number is around four as expected. The derivative of the coordination number is more illustrative [fig¬ 
ure [3] (b)]. It shows the location of the first, the second shell and even the third one, the latter being 
developed at very low temperatures and serving as a manifestation of a spatial order that develops in the 
system and involves a larger number of particles than the first neighbors. 

To proceed, it is necessary to introduce the concept of H-bonds. Energetic and geometric criteria, un¬ 
avoidably involving certain degree of arbitrariness, are used in the hterature to interpret the structure 
of water and aqueous solutions. A comprehensive analysis of several geometric definitions of hydrogen 
bonding, leading to different values for the average number of H-bonds per molecule, has been per¬ 
formed in reference (sj]. 

In this work, we use a single geometric criterion, see e.g. Zhang et al. jsil] for two different water mod¬ 
els. Three conditions determine the existence of the H-bond, namely the distance Rqq between oxygen 
atoms belonging to two water molecules should not exceed the threshold value the distance i? 0 H 
between the acceptor oxygen atom and the hydrogen atom connected to the donor oxygen should not 
exceed the threshold value R^^. Finally, the angle 0-H-- 0 should be smaller than the threshold value, 
see e.g. references concerning the analysis of the structure of liquid methanol as well as of su¬ 

percritical water. If the coordinates of two molecules fulfill the above conditions, they are considered as 
H-bonded. Zhang et al. proposed constant threshold values for the distances, 3.S A and 2.4S A, for R^^ 
and respectively. However, the values used in this work are the ones coming from the correspond¬ 
ing radial distribution functions minima for 0-0 and 0-H at each calculated temperature. In contrast to 
Jszll . we used the threshold constant value for the angle 6, complementary to the angle a between the 
vectors OH and rou (see figure 1 of reference (^), that is equal to 150°. With this definition, we obtained 
the average number of hydrogen bonds per water molecule similar to the value by Bako et al. who 
used the energetic criterion for bonding. 

The temperature dependence of the number of H-bonds averaged upon simulation time and normal¬ 
ized by the number of water molecules in the system, {«hb)jVi„. following from our calculations for two 
water models, the SPC/E and TIP4P-Ew, is shown in figure |4] (a). This figure explains that a decreasing 
temperature results in a larger average number of H-bonds. It can be seen that there is a change of the 
slope at r =: 180 K for the SPC/E model where the transition into a glassy state is reported. On the other 
hand, the change of the slope is observed at a higher temperature, around T - 200 K, for the TIP4P-Ew 
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bonding state of a water molecule T(K) 


Figure 4. (Color online) Temperature dependence of the average number of H-bonds per water molecule 
(panel a). Distribution function of tbe angles formed by three oxygens on different molecules at different 
temperatures (panel b). Fractions of water molecules in different bonding states on temperature (panel 
c). A comparison of fractions of differently bonded molecules on temperature for the SPC/E and TIP4P-Ew 
models (panel d). 


model. However, this latter model is characterized by a higher number of H-bonds per molecule in the en¬ 
tire temperature interval studied. In both cases we see that the average number of bonds does not reach 
four, or, in other words, even at low temperatures, perfect tetrahedral configuration is not attained. 

The same conclusion can be reached by considering the “bond-angle distribution function” already 
studied by Bryk and Haymet (3- It describes the probability to find certain configurations involving 
three oxygens belonging to different water molecules. The threshold distance is assumed to be at the first 
minimum of the pair distribution funcion, goo (f) at each temperature studied. The data were normalized 
by the number of molecules that one finds in the shell limited by the threshold distance. All the curves 
given in figure|4](b) are characterized by two maxima. One of them at 0 =; 53° describes the configurations 
characteristic of a simple liquid (3- On the other hand, the second maximum at 6 =; 107° is close to the 
angle 0 =; 109.5° describing a perfect tetrahedral configuration. At a high temperature, T = 360 K, a simple 
liquid type structure prevails. While the system is cooled, the first maximum decreases and becomes very 
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small at r = 200 K, whereas the maximum describing the configurations close to tetrahedral substantially 
grows in value. It is worth mentioning that both maxima slightly shift to higher angles with a decreasing 
temperature. Still, even at low temperatures, the shape of the curves is not delta-like, witnessing a certain 
degree of disorder in the system. Moreover, we do not observe any peculiarities of this property either at 
TMD or at the temperature where the density is at minimum. 

In order to obtain a better insight into the changes of the structure in terms of bonds with temper¬ 
ature, we calculate the time-ave rag es of fractions of differently bonded molecules. This representation 
was already used by Bako et al. ISSh in their analysis of the simulated structures of water-methanol mix¬ 
tures at room temperature. Our results concerning the SPC/E model at different temperatures are shown 
in figure|4](c). We observe that at a high temperature, T - 330 K, the fractions of molecules participating 
in two and three bonds are the largest. Comparing this and room temperature, it can be seen that the 
fraction of molecules with four bonds grows at the expense of a decreasing fraction of molecules with 
a single bond and two bonds. The fraction of particles with three bonds remains practically unchanged. 
A further decrease of temperature down to T = 240 K (note that it is very close to the TMD) leads to a 
substantial increase of four-bonded particles again at the expense of singly-bonded and doubly-bonded 
species. The fraction of triple-bonded species remains intact. If the temperature falls down to 200 K, not 
only the fraction of 4-bonded particles grows substantially, but the fraction of triple-bonded species ex¬ 
hibits a strong decrease in value. To summarize, the growth of density is accompanied by the growth of 
the fraction of 4-bonded particles with almost constant fraction of triple-bonded species. On the other 
hand, the descending density region, below the TMD, is characterized by the growth of 4-bonded species 
with a simultaneous decrease of fractions of triple-bonded, double-bonded and single bonded-molecules. 
Interestingly, the changes of the structure at even lower temperature are very small (we do not show 
these curves to avoid overload of the figure). 

The dependence of fractions of differently bonded molecules on temperature, in the entire tempera¬ 
ture range studied, is shown in figure|4](d). Here, we present our results for two different models, namely 
for the SPC/E and TIP4P-Ew. The pattern is very similar in both cases. The TMD for each model obtained 
from the calculations of density on temperature is located very close to the crossing points between 
the curves describing the dependence of triple-bonded and 4-bonded molecules. The TMD values are 
marked in the figure as well. It is worth mentioning that the curves describing the fraction of triple- 
bonded molecules is not sensitive to the temperature variable in the region between high temperatures 
and T =; 240 K. Moreover, the values for this fraction coming from the calculations of two models in ques¬ 
tion practically coincide in an ample interval of temperatures. However, the difference between these 
two models is well observed due to their capability of forming structures describing 4-bonded species. 

In order to complete our discussion of the results given in all panels of figureH) we would like to make 
the following comments. As mentioned above, the geometric definition of H-bonds permits flexibility in 
the sense that one can attempt to change the values and the set of parameters of H-bond definition. In this 
respect, the energetic definition is more restrictive. Besides the choice of a set (ifoo. ^'oh. d), we also tested 
the sets (Rqo, ron) and (ifoo. P, cf. reference (sj) figure 1). These less restrictive geometric definitions lead 
to the results quahtatively similar to the ones discussed above. Namely, the fraction of quadruply-bonded 
species grows with a decreasing temperature whereas all other fractions decrease. Unfortunately, we 
have not observed any characteristic crossing attributable to the TMD with these definitions. Therefore, 
we hope to extend our study of the models in question in the framework of energetic definition of H- 
bonding as weU. 

The properties described above refer to the structure of the systems in question mostly at short inter¬ 
particle separation. On the other hand, the long-range, asympthotic behavior of correlations between 
molecules possessing dipole moment is described by the dielectric constant, e, see e. g. 13911 . In general 
terms, the calculation of the dielectric constant from simulations is a difficult task (2i,l4(il4lll . especially 
at low temperatures. Several factors can influence the result, such as e.g., the number of molecules, the 
type of thermostat and barostat, precision of the long-range interactions summation. Actually, very long 
runs are necessary to obtain reasonable estimates for this property, as it is usually calculated from the 
time-average of the fluctuations of the total dipole moment of the system as follows: 

( 1 ) 

ShTV ^ ' I’ 
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Figure 5. (Color online) Time dependence of the dielectric constant of the SPC/E and TIP4P-Ew water 
model at different temperatures from our NPT molecular dynamics simulations, panels a and b, respec¬ 
tively. Time dependence of the Cartesian projections of the total dipole moment of the system at room 
temperature for the SPC/E model (panel c). Temperature dependence of the dielectric constant for the 
SPC/E and TIP4P-Ew models from our calculations, from the fitting procedure by Fennell et al. and 
the experimental data 


where is the Boltzmann constant and V is the simulation cell volume. Some of our results at different 
temperatures for the SPC/E model and for TlP4P-Ew are shown in figure[3(a) and[3(b), respectively. We 
would like to recall here that all the data given in figure were obtained for the systems having two 
thousand particles. Practically all the runs lasted more than 20 ns. 

Still, as it can be seen, there is room for experimenting in order to get better results below the room 
temperature by extending the length of the runs or considering larger systems. In addition, we plan to ex¬ 
plore the effects of changing the cutoff distance, of including the long-range corrections and the precision 
of Ewald summation as well as to test other thermostats. All the results for each model, however, refer to 
the temperatures above the respective melting point. Still, there is one delimiting factor. It manifests itself 
in the polarization of the samples at low temperatures. At room temperature, figureO(c), the projections 
of the total dipole moment of the system tend to zero at large times, whereas at lower temperatures, in 
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the supercooled regime and close to the formation of glassy states, the samples inevitably fall into con¬ 
figuration with non-zero one or more Cartesian projections of M. Then, the calculations of the dielectric 
constant require different algorithms. A summarizing insight into our data for e{T) for two models in 
question are given in figure[^(d). In general terms, the inclination of the temperature dependence of the 
dielectric constant is correctly described by two models, compared to the experimental results. The fitting 
of the dielectric constant at room temperature does not garantee the correct temperature dependence, as 
we can see from the results of the modified model jilll . 

To summarize this brief report, we have used molecular dynamic simulations in the NPT ensemble 
to study the density anomaly of water for the SPC/E and TlP4P-Ew models and to establish the relation 
between this anomaly with the microscopic structure of the system. The structure is described in terms of 
radial distribution functions, oxygen-oxygen coordination number and its derivative, the average num¬ 
ber of H-bonds per molecule, angular distributions and fractions of molecules characterized by a differ¬ 
ent number of bonds. We observed that at temperature of maximum density, the fractions of molecules 
with three and with four bonds are almost the same, or in other words it is a reasonably good structural 
estimator for the TMD. In addition, we presented our preliminary results concerning the temperature 
dependence of the dielectric constant and compared them with experimental data. While the slope of 
this dependence is weU described, the values are underestimated compared to the experimental results. 
It would be desirable in a future work to extend our structural analysis to find estimates for the cooper¬ 
ative bonding of molecules and relate them with the density and other anomalies of water. On the other 
hand, it would be interesting to estabhsh the relationship between the structural properties describing 
the electric properties of the system and the dielectric constant. 
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TeMneparypHa 3a,ne>KHicTb MiKpocKoniMHoV CTpyKTypi/i i 
aHOMa/iifl rycTi/iHi/i b MOAe,n5ix boai/i SPC/E ra TIP4P-Ew. 
Pe3y/ibTaTi/i m 0 ^ 6/1 lOBaHHH mbtoaom MO/ieKy/inpHOi Ai^HaiviiKi/i 

E. raBicia-AHflpecPI f. /\0MiHreP, 0. fliBicP 

^ iHCTi/nyT xiMii, HaqiOHanbHuii aBTOHOMHi/it) yHiBepcuTer Mbkcmkh, MexiKO, MexcuKa 
^ iHCTi/nyT flOcniflxeHHfl MarepianiB, HaqiOHanbHuii aBTOHOMHuii yHiBepcMrer Mbkcmkh, Mexixo, MexcuKa 

Mi/i floc/iifli/i/m reMnepatypHi aa/iexHOcri MiKpocKoniuHOi crpyKTypn MOfle/ieCi BOfli/i SPC/E i TIP4P-Ew b rep- 
Minax napHMX c|)yHKL(iii poanofliny, KOopflUHaqiiiHnx nucen, cepe^Hboro nucna BOflHeBnx SB'naKiB, poanofliny 
3B'naaHnx CTaniB OflHiei MoneKy/m, BHKopncTOByK)HH Merofl MO/ieKynnpHOi flunaMiKH npw nocritiHorMy Ti/icxy. 
EBO/iKDqin crpyKTypn nocraB/ieHa y BiflnoBiflHiCTb flo aanexHOcri rycn/iHn BOfli/i Bifl reMneparypn b o6nacTi, tqo 
3HaxoflnTbcn Mix bucokhmi/i reMneparypaMn i TeMneparypaMi/i, fle cncreMa crae nepeoxo/io/pKeHOKD. flOKaaa- 
HO, tpo uacTKa MO/ieKy/i a rpbOMa i nornpMa aB'^axaMi/i Bnauauae Maxci/iManbHy rycn/iHy flnn o6ox MOflenen. 
Binbtu Toro, orpi/iMaHO i npoaha/iiaoBaHO reMneparypHy aa/iexHicrb flie/ieKipi/iuHOi cranoi. 

KnioMOBi cnoBa: MOflejii BOfii/i, aHOMajiia rycri^Hi^, MoneKynapna fiMnaMiKa 
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